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AN ALGORITHM FOR INTEGRATING 
SIMULTANEOUSLY THE RESTRICTED PROBLEM AND 
ITS FIRST VARIATIONS IN THIELE'S COORDINATES 

R. H. Estes 


ABSTRACT 


Iii the search for periodic orbits in the planar restricted prob- 
lem of three-bodies, it is of interest to integrate simultaneously 
the orbital and variational equations in regularizing coordi- 
nates . 2 Estes and Lancaster 5 have presented an efficient algo- 
rithm for the integration of the equations of motion of the 
restricted problem in Thiele's variables and it is the purpose 
of this note to extend the algorithm to include the variational 
equations. 
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AN ALGORITHM FOR INTEGRATING 
SIMULTANEOUSLY THE RESTRICTED PROBLEM AND 
ITS FIRST VARIATIONS IN THIELE'S COORDINATES 


EQUATIONS OF MOTION AND VARIATION 

The Lagrangian for the plar xr restricted three-body problem in Thiele coordi- 
nates with the origin at the midpoint between the primaries and the x-axis passing 
through both primaries while rotating with the mean motion may be written 

S. ~ j(u' 1 2 + v‘ 2 ) - g(u' sinh 2v + v' sin 2u) D (1) 


where 


fi 


= ~ (cosh v - [2 fj.- l] cos u ) + ~ U j cos u cosh v - J - j [cosh 2 v + cos 2 u]j 


t' - ^(cosh 2 v- cos 2 u) - fjfj 


1 

fj - j ( c °sh v - cos u) 


1 

r 2 = ^ (cosh v + cos u) 


u' 2 + v' 2 - 2 1) 


r 1 r 2 


(1-, )r{ 


2 + M r 2 


2 ( 1 - M ) + 2jj. 
r i r 2 


and where y is the ratio of the mass of the smaller primary (assumed to be to 
the left of the origin) to the sum of the two primary masses. Here prime ac- 
cents indicate differentiation with respect to a regularizing pseudo-time r, r t 
and r 2 represent the distances between the primaries and the infinitesimal third 
body and J represents the Jacobi constant of motion. The variables u and v 
are related to the Cartesian coordinates by the transformation 


1 



x - 


2 cos u cosh 


y - ~ ~2 sin u sinh v 


The Lagrangian equations of motion become 


( 2 ) 


u" - 2t'v* " 


v" + 2t*u’ = 


bo 

Bu 


bh 

2 v 


( 3 ) 


where 


3 D sin u 

3 u 16 


[8(2/i - 1) - 4J cos u + (2/x - 1) cosh v(3cos 2 u - cosh 2 v) + 2cos 3 u] 


3 n sinh v r , - 

yy - jg — [8 * 4J cosh v *■ (2/x - 1 ) cos u (3cosh* v - cos 2 u) + 2 cosh 3 v J . 


The solution of Equations (3) depends upon initial conditions and the parameters 
A* and J. Thus 


U = u(r; u 0 , V 0 , u^, Vj, fX, J) 
v = v(t; u 0 , V 0 , V-, fi, J) . 


3 u d v 

lenoting u e = yy, v e = yy, etc,, where e represents one of the parameters or 
nitial conditions, the first variational equations of the system (3) become 

!u^ - 4 1 ’ vj = u e | cos u [(2 ix ~ 1) “ y cos u + (2/x - 1) cosh v(3cos 2 u ~ cosh 2 v) 
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+ 2 cos 3 u + 2 sin u v' + — cosh v (cos 2 u - 1) 


3 , 

+ COS 3 


-l-MH 3 ') 


cos 2 u 


v e |sinh v ^2v' cosh v - — ™ — — sin u (cosh 2 v - cos 2 u)j| 


sin u cosh v , , 

+ ix e ' 2 sinu + ^ (3 cos 3 u - cosh 3 


•>} 


~2 J e sin u cos u 


2v" + 4t' u; = v e 


K- 


cosh v 1 - ~ cosh v + (2/i - 1) cos u (3 cosh 2 v - cos 2 u) 


, 3 2u “1) 

+ 2 cosh 3 v - 2 sinh v u' + ■ - — ^ cos u (cosh 2 v - 1) 

+ -|cosh 3 vj + yj - ( 3 + 4 2 ^ ) cosh 2 v| 

+ u e jsin u j^ - 2 u‘ cos u ^ sinh v (cosh 2 v - cos 2 u) | 


{ sin u sinh v 
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(3 cosh 2 v - cos 


«>1 


A 

2 J g sinh v cosh v 


(4) 


POWER SERIES SOLUTION 


We now extend the notation introduced in Reference [4] . Let 



J K 1 

k, = 2fj. - 1 , k 2 = - - 7 *. k 3 = - Y' k 4 = ~ 3k 3 (5) 

a = sinu, b = cos u, c - sinh v, d - cosh v ( 6 ) 

p “ b 2 , q = d 2 , w = q - p (7) 

r = k 4 d + b/4, s = k 4 b + d/4 ( 8 ) 

f = kj + k 2 b + rp + k 3 dq (9) 

g = l + k 2 d+$q + k 3 bp (10) 

u' = a, v 1 = / 3 , Uj = a e , v' t ~ (11) 

y = f + 2 a /3 * — bp + 2 k 4 d(p - 1 ) ( 12 ) 

\ = 2y3d-k 4 aw (13) 

p - g-2ca+-|-dq+ 2k 4 b(q - 1) (14) 

a = - 2 b a - k 4 cw (15) 

h = d(3p - q), z = b(3q - p) (16) 

r = " k 2 “ (| ~ k 2 ) P + by (17) 

A = c\, S = acr (18) 

R = " k 2 " (I " k 2 j q + dp (19) 

H = a (2 + , L = | z (20) 
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Equations (3) and (4) then become 

2 a' = af+w/3 (21) 

2/3* = cg-wa (22) 

2a' e ~ T + v e A + w + p e H - j J e a b (23) 

2 fig = v e R + u { S“wa ( cd, (24) 


Assume the solutions of (7) - (24) can be represented by the power series 


u 


£ u ‘ t1 ' 
i = 0 


V 



co 



i *0 



i =0 i =0 


and similar series for a, b, c, d, p, q, w, r, s, f, g, y, p, a-, h, z, P, A, 
R , S, H, and Z in a neighborhood of r = 0. Substituting these series into 
(7) - (24) and equating the coefficients of each power of r gives, for i > 0 



j*0 j*0 


r i = k 4 d i f s i - k « b i t- dj/4 (26) 


f. 


k l S i0 * k 2 b i 



Pi-i 


+ k, 



d q _ . 


(27) 
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(28) 


«i = 8 

2(i 
2(i 

(i + 2 > u i + 2 
?i = f » - 2k 4 d t 

Pi = 8i " 2k 4 b i ' 

N - 

j-o 

i 

h i = ^ 2] d i Pi-j 

j *o 

r. = 

i 

R. = 

i 


iO + k 2 d i 




> 

+ 1) &i+i = Z] <8j fi ‘ j + Wj A-i> 

j-0 


1 

l > = Z] (Cj gi * J ’ W i ai -J } 

j-0 


(29) 


(30/ 


i + i* 


(i + 2) v i + 2 = /3. + 1 , 4(i + i) * » w. (31) 



*J Pi'i 



b i Pi-) * 2k, d ) Pl-i < 32 > 

j-0 






b. 


j *»»-j 


(33) 





Wj c i . j (34) 



-Mio- (t* k i) Pi + 22 b|7 ‘- j 

j -0 

-Mio - ♦ 22 d / p i-i 

j-0 


(36) 


(37) 
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(38) 


i i 

L i = 2L1 Si S Zj a » 

j-0 j-0 


H = 2 a 


1 i i 

i + 4 y] ®i h l-j' Z i 4 C i *i“j 

j-0 j-0 


(39) 


1 i 

j-0 j-0 


(40) 


1 1 

= ' i « Z i*2Z [ v «, IR l-i *<Vr-i * u <, 3 l-il - 2 J < X! C > d -i (41) 

j-0 


i*0 


(i + 2) u = a . (i + 2') v = /3 

v * «!♦! «i*r v *' v «i*a ^«i+t 


(42) 


X X 

< i + 1 > a iM = 7] bj ai "j’ (i + 1)bi + 1 = ~T] a > a '-’ 

j-0 j-0 


( 43 ) 


( i + i > c i.i = 2Z (i + i)d **' = l c - *-< (44) 

j-0 j-0 


Here 5. 0 Is the Kronecker symbol which equals one if i is zero, and zero 
otherwise. Given the initial conditions u Qt v 0 , a Q - u t , 0 O = v t , t 0 , u e , v t , 
a * u <t and “ v e ,» (0) may be solved for a 0 , b 0 , c 0 , d 0 and the algo- 
ritnm (25) - (44) will then give all higher order coefficients in terms of the pre- 
ceding coefficients, so that the r lution may be extended optimally by analytic 
continuation using a variable step size and a variable number of terms in the 
power series expansion for each integration step. 4 It is to be noted that the 
equations of motion may be computed independently by deleting Equations (32) - 
(42) from the above algorithm. 
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As an example of applic ability this algorithm may be used In conjunction with 
the well-known predictor-corrector technique for the numerical computation of 
natural families of periodic orbits. 1,2 Following this method we set = 0 and 
the predictor variational equations would have J, =1, resulting in the expected 
nonhomogeneous system of linear equations. The iterative correcting procedure 
is isoenergetic so that J e = 0 and the homogeneous system of variational equa- 
tions yielding the corrected periodic orbit gives the characteristic stability 
exponents as eigenvalues of the numerically integrated matrizant at the end of 
the fundamental period. 3 
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